Bose-Einstein condensation of quasiparticles in graphene 
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■ The collective properties of different quasiparticles in various graphene based structures in high 

04 ' magnetic field have been studied. We predict Bose-Einstein condensation (BEC) and superfluidity 

of 2D spatially indirect magnetoexcitons in two-layer graphene. The superfiuid density and the 
!—{ , temperature of the Kosterlitz-Thouless phase transition are shown to be increasing functions of 

the excitonic density but decreasing functions of magnetic field and the interlayer separation. The 
instability of the ground state of the interacting 2D indirect magnetoexcitons in a slab of superlattice 
with alternating electron and hole graphene layers (GLs) is established. The stable system of 
04 ' indirect 2D magnetobiexcitons, consisting of pair of indirect excitons with opposite dipole moments, 

is considered in graphene superlattice. The superfiuid density and the temperature of the Kosterlitz- 
Thouless phase transition for magnetobiexcitons in graphene superlattice are obtained. Besides, the 
BEC of excitonic polaritons in GL embedded in a semiconductor microcavity in high magnetic field is 
predicted. While superfiuid phase in this magnetoexciton polariton system is absent due to vanishing 
of magnetoexciton-magnetoexciton interaction in a single layer in the limit of high magnetic field, 
. the critical temperature of BEC formation is calculated. The essential property of magnetoexcitonic 

systems based on graphene (in contrast, e.g., to a quantum well) is stronger infiuence of magnetic 
field and weaker infiuence of disorder. Observation of the BEC and superfiuidity of 2D quasiparticles 
in graphene in high magnetic field would be interesting confirmation of the phenomena we have 
' described. 

PACS numbers: 03.75.Hh, 73.20.Mf, 71.36.-fc 

O ' I- INTRODUCTION 

The production of graphene, a two-dimensional (2D) honeycomb lattice of carbon atoms that form the basic planar 
structure in graphite, has been achieved recently [l|,lj| . Electronic properties of graphene, caused by unusualproperties 
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K*" ' of the band structure, became the object of many recent experimental and theoretical studies [l|, 
^} Graphene is a gapless semiconductor with massless electrons and holes, described as Dirac-fermions [8L The various 
studies of unique electronic properties of graphene in a magnetic field have been performed recently fiol fill . fT3 | . 
: The energy spectrum and the wavefunctions of magnetoexcitons, or electron-hole pairs in a magnetic field, in graphene 
have been calculated in interesting works [H, • 

The 2D electron system was studied in quantum wells (QWs) [l^. The systems of spatially- indirect excitons (or 
pairs of electrons and holes spatially separated in different QWs) in th e sy stem of coupled quantum wells (CQWs), 
with and without a magnetic field were studied in Refs. [16|, |17|, [11 m US mi, m, m. The experimental and 
theoretical interest to study these systems is particularly caused by the possibility of the BEC and superfluidity of 
. 5^ ^ indirect excitons, which can manifest in the CQWs as persistent electrical currents in each well and also through 
^ . coherent optical properties and Josephson phenomena [la . [Tsl . [l9l [20I [2^ . The outstanding experimental success 



' was achieved now in this fie ld [2^, [25l . [26l. \2l\. The electron-hole pair condensation in the graphene-based bilayers 
5^ , have been studied in [H, [H, l30l. l3l| . 

The collective properties of Bose quasiparticles such as excitons, biexcitons, and polaritons in various graphene- 
based structures in high magnetic field are very interesting in the relevance to the BEC and superfluidity, since the 
random field in graphene is weaker than in a QW, particularly, because in a QW the random field is generated due 
to the fluctuations of the width of a QW. Let us mention that if the interaction of bosons with the random field is 
stronger, the BEC critical temperature is lower [s^ . In this paper we study the superfluidity of magnetoexcitons 
in bilayer graphene, instability of the system of magnetoexcitons in superlattice formed by many GLs, superfluidity 
of magnetobiexcitons in graphene superlattice, and BEC of polaritons in GL embedded in optical microcavity in a 
trap. Let us mention that all these systems of quasiparticles are considered in high magnetic field. The BEC of 
magnetoexcitons in graphene structures can exist at much lower magnetic field than in QWs, because the distance 
between electron Landau levels in graphene is much higher than in a QW at the same magnetic field, and, therefore, 
the lower magnetic field is required in graphene than in a QW to neglect the electron transitions between the Landau 
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levels. 



II. EFFECTIVE HAMILTONIAN OF MAGNETOEXCITONS AND PHOTONS IN MICROCAVITY IN 

HIGH MAGNETIC FIELD 



Recently, Bose coherent effects of 2D exciton polaritons in a quantum well embedded in an optical microcavity 
have been the subject of theoretical [ssj . and experimental (3^ . [ssl . [36| studies. To obtain polaritons, two mirrors 
placed opposite each other form a microcavity, and quantum wells are embedded within the cavity at the antinodes 
of the confined optical mode. The resonant exciton-photon interaction results in the Rabi splitting of the excitation 
spectrum. Two polariton branches appear in the spectrum due to the resonant exciton-photon coupling. The lower 
polariton (LP) branch of the spectrum has a minimum at zero momentum. The effective mass of the lower polariton is 
extremely small, and lies in the range 10~^ — 10~^ of the free electron mass. These lower polaritons form a 2D weakly 
interacting Bose gas. The extremely light mass of these bosonic quasiparticles, which corresponds to experimentally 
achievable excitonic densities, result in a relatively high critical temperature for superfluidity, of 100 K or even higher. 
The reason for such a high critical temperature is that the 2D thermal de Broglie wavelength is inversely proportional 
to the mass of the quasiparticle. 

While at finite temperatures there is no true BEC in any infinite untrapped 2D system, a true 2D BEC can exist 
in the presence of a confining potential (stI . [ssj . Recently, the polaritons in a harmonic potential trap have been 
studied experimentally in a GaAs/AlAs quantum well embedded in a GaAs/AlGaAs microcavity [s^. In this trap, 
the exciton energy is shifted using stress. In this system, evidence for the BEC of polaritons in a QW has been 
observed [l^l . The theory of the BEC and superfluidity of excitonic polaritons in a QW without magnetic field in 
a parabolic trap has been developed in Ref. |41ii|. The Bose condensation of polaritons is caused by their bosonic 
character [iol [4ll |4^ . However, while the exciton polaritons have been studied in a QW, the formation of the 
polaritons in graphene in high magnetic field have not yet been considered. Moreover, the magnetopolaritons formed 
as superposition of magnetoexcitons and cavity photons in magnetic field have not yet been studied. We consider a 
2D system of polaritons in graphene layers (GLs) embedded in a microcavity from the point of view of the existence 
of the BEC within it. 

Lets us consider the most general case when the superlattice with alternating electronic and hole parallel GLs in 
the external field is embedded in an optical microcavity in high magnetic field. At the small densities n the system 
of indirect excitons at low temperatures is the two-dimensional weakly nonideal Bose gas with normal to wells dipole 
moments d in the ground state (d = eZ3, e is the charge of an electron, D is the interlayer separation). In contrast to 
ordinary excitons, for the low-density spatially indirect magnetoexciton system the main contribution to the energy is 
originated from the dipole-dipole interactions C/_ and C/+ of magnetoexcitons with opposite (see Fig.[T|) and parallel 
dipoles, respectively. The potential energy of interaction between two indirect magnetoexcitons with parallel U+{R) 
and opposite U- (R) dipoles is a function of the distance R between indirect magnetoexcitons along GLs and is given 
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2p2 p2 

U-{R) = + , (1) 

where e is the charge of an electron, e is the dielectric constant. 

The Hamiltonian of magnetoexcitons and photons in the strong magnetic field is given by 

Hfot — ^mex ^ -^ph ^" ^mex—ph ; (2) 

where Hmex is a magnetoexcitonic Hamiltonian, Hph is a photonic Hamiltonian, and H^xc-ph is the Hamiltonian of 
magnetoexciton-photon interaction. 

Let us analyze each term of the Hamiltonian The effective Hamiltonian and the energy dispersion for magne- 
toexcitons in graphene layers in a high magnetic field B in the infinite system was derived in Ref. [43]. The effective 
Hamiltonian Hmex of the low-density system of the indirect magnetoexcitons in high magnetic field in the superlat- 
tice in the subspace of the lowest Landau level is given by [isj (we neglect the electrons transitions between different 
Landau levels due to electron-hole Coulomb attraction) 

Hmex ^ Hq + Hint ■ (3) 
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FIG. 1: Two-dimensional indirect magnetobiexcitons consisting of indirect magnetoexcitons with opposite dipole moments, 
located in neighboring pairs of GLs. 

Here Hq is the effective Hamiltonian of the system of noninteracting trapped magnetoexcitons in high magnetic field: 

-H'o = ^ ernex{p){apap + b^bp + atpa_p + btpb_p) , Smex{P) = Eband{r) - + £o{P) , (4) 

p 

where tip , dp , 6p , b-p 8,rc creation and annihilation operators of the magnetoexcitons with up and down dipoles. In 
Eq. Ebandif) is the band gap energy which can depend on the position of the magnetoexciton in space in the 
presence of the trap, is the binding energy of a magnetoexciton, and £q{P) = P^ /{2mB), where mg is the effective 
magnetic mass of a magnetoexciton. Similarly to the Bose atoms in a trap in the case of a slowly varying external 
potential [i^l, we can make the quasiclassical approximation, assuming that the effective magnetoexciton mass does 
not depend on a characteristic size / of the trap and it is a constant within the trap. This quasiclassical approximation 
is valid if P 3> h/l. The harmonic trap is formed by the two-dimensional planar potential in the plane of graphene. 
The potential trap can be produced in two different ways. One way is when the potential trap can be produced 
by applying an external inhomogencous electric field. The spatial dependence of the external field potential V{r) 
is caused by shifting of magnetoexciton energy by applying an external inhomogencous electric field. The photonic 
states in the cavity are assumed to be unaffected by this electric field. In this case the band energy is given by 
Eband{r) = Eband{0) + V'(r) {Eband{0) = \/2hv F / r B is the band gap energy, which is the difference between the 
Landau levels 1 and 0, vp is the Fermi velocity of electrons in graphene [1^, tb — \/hc/JeB) is the magnetic length). 
Near the minimum of the magnetoexciton energy, V{r) can be approximated by the planar harmonic potential 7r^/2, 
where 7 is the spring constant, r is the distance between the center of mass of the magnetoexciton and the center of 
the trap. Note that a high magnetic field does not change the trapping potential in the effective Hamiltonian [i^. Let 
us mention that, while the quasiparticles in GLs and QWs in high magnetic field are described by the same effective 
Hamiltonian with the only difference in the effective magnetic mass of magnetoexciton. This difference is caused by 
the four-component spinor structure of magnetoexciton wave function in GL, while magnetoexciton in a QW and 
CQWs is characterized by the one-component scalar wave function. 

The Hamiltonian which describes the interaction between magnetoexcitons is 

H^nt = (2S')-i Epi-HP2=P3+P4(f^+(ap,ap3aP2api + b+b+bp^bp^) - U ^a+b+^ap^bp^), (5) 

where (7+ and U- are the 2D Fourier images of U+{R) and U-{R), respectively, and S is the surface of the system. 
The Hamiltonian and the energy spectrum of non-interacting photons in a semiconductor microcavity are given 



4 



by (47|: 



where Op and ap are photonic creation and annihilation Bose operators. In Eq. Lq is the length of the cavity, 
n = y/ec is the effective refractive index and ec is the dielectric constant of the cavity. We assume that the length of 
the microcavity has the following form: 

Lc{B) = — ^ , (7) 

corresponding to the resonance of the photonic and magnetoexcitonic branches at P = 0, i.e. emex{0) — £ph{0)- 
As it follows from the energy spectra in ^ and and Eqs. ^ and ([7]), the length of the microcavity, corre- 
sponding to a magnetoexciton-photon resonance, decreases with the increment of the mag netic field as B-^/^. The 
resonance between magnetoexcitons and cavity photonic modes can be achieved either by controlling the spectrum of 
magnetoexcitons s^xiP) by changing magnetic field B or by choosing the appropriate length of the microcavity Lq- 
Alternatively to the case considered above, now the trapping of magnetopolaritons is caused by the inhomogeneous 
shape of the cavity when the length of the cavity is determined by Eq. ([7]) with the term added to i?band(0)— fjj^'', 

where r is the distance between the photon and the center of the trap. In case 2, the 7 is the curvature characterizing 
the shape of the cavity. The Hamiltonian and the energy spectrum of the photons in this case are shown by Eq. ©, 
and the length of the microcavity is given by Eq. ([7]). In this case, for the photonic spectrum in the effective mass 
approximation is given by substituting the slowly changing shape of the length of cavity depending on the term 7r^/2 
into Eq. ^ representing the spectrum of the cavity photons. This quasiclassical approximation is valid if P H/l, 

where I = {h/ {uib'^o))^^'^ is the size of the magnetoexciton cloud in an ideal magnetoexciton gas and luq — \plj™^- 
The Hamiltonian and energy spectrum of magnetoexcitons in this case are given bv (|4]). 

The Hamiltonian of the magnetoexciton-photon coupling has the form (see Refs. (48l l49l [sol [5l|): 

Hmex-ph = hQ.R ^ ap6p + h.c. , (8) 
p 

The projection of the electron-hole Hamiltonian in magnetic field for the CQWs and GLs onto the lowest Landau 
level results in the effective Hamiltonian (|15p with renormalized mass and where term related to the vector potential 
is missing. The magnetic field in the effective Hamiltonian (|15p enters in the renormalized mass of the magnetoexciton 
niB- Therefore, Hamiltonian for the spatially separated electrons and holes in two-layer system for the CQWs and 
for the bilayer graphene can be reduced in high magnetic field to the effective Hamiltonian (fT5|) . Magnetic field B is 
reflected by the effective Hamiltonian (|15p only through the effective magnetic mass of a magnetoexciton ms in the 
expression for £o(P) in the first term of Hmex- The only difference in the effective Hamiltonian (fT5| for the CQWs 
and bilayer graphene realizations of two-layer systems is that for the bilayer graphene is four times less than for 
the CQWs due to the four-component spinor structure of the wavefunction of the relative motion for the isolated 
non-interacting electron-hole pair in magnetic field [52l| . 

Transitions between Landau levels due to the Coulomb electron-hole attraction for the large electron-hole separation 
D ^ rs can be neglected, if the following condition is vahd: Eb — e^/(eb-D) <^ hwc = heB{me + mh) I (2memhc) 



for the QWs and Ei, = Ae^/^eD) <C hvp/rB for the GLs, where vp is the Fermi velocity of electrons [45|, Ei, and 
LUc are the magnetoexcitonic binding energy and the cyclotron frequency, respectively. This corresponds to the high 
magnetic field P, the large interlayer separation D and large dielectric constant of the insulator layer between the 
GLs. 



III. BOSE-EINSTEIN CONDENSATION OF POLARITONS IN GRAPHENE IN A HIGH MAGNETIC 

FIELD 

In this Section we consider trapped polaritons in a single graphene layer embedded into an optical microcavity in 
high magnetic field. When an undoped electron system in graphene in a magnetic field without an external electric 
field is in the ground state, half of the zeroth Landau level is filled with electrons, all Landau levels above the zeroth 
one are empty, and all levels below the zeroth one are filled with electrons. We suggest using the gate voltage to control 
the chemical potential in graphene by two ways: to shift it above the zeroth level so that it is between the zeroth 
and first Landau levels (the first case) or to shift the chemical potential below the zeroth level so that it is between 
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the first negative and zeroth Landau levels (the second case). In both cases, all Landau levels below the chemical 
potential are completely filled and all Landau levels above the chemical potential are completely empty. Based on 
the selection rules for optical transitions between the Landau levels in single-layer graphene [53|, in the first case, 
there are allowed transitions between the zeroth and the first Landau levels, while in the second case there are allowed 
transitions between the first negative and zeroth Landau levels. Correspondingly, we consider magnetoexcitons formed 
in graphene by the electron on the first Landau level and the hole on the zeroth Landau level (the first case) or the 
electron on the zeroth Landau level and the hole on the Landau level —1 (the second case). Note that by appropriate 
gate potential we can also use any other neighboring Landau levels n and n + 1. 

For the relatively high dielectric constant of the microcavity, e ^ /(hvp) ~ 2 the magnetoexciton energy in 
graphene can be calculated by applying perturbation theory with respect to the strength of the Coulomb electron- 
hole attraction analogously as it was done in [iTj for 2D quantum wells in a high magnetic field with non-zero electron 
•nie and hole m/j masses. This approach allows us to obtain the spectrum of an isolated magnetoexciton with the 
electron on the Landau level 1 and the hole on the Landau level in a single graphene layer, and it will be exactly 
the same as for the magnetoexciton with the electron on the Landau level and the hole on the Landau level — 1. 
The characteristic Coulomb electron- hole attraction for the single graphene layer is e^/{erB)- The energy difference 
between the first and zeroth Landau levels in graphene is hvp/rs- For graphene, the perturbative approach with 
respect to the strength of the Coulomb electron-hole attraction is valid when /{ers) ^ hvp/rs [13 • This condition 
can be fulfilled at all magnetic fields B if the dielectric constant of the surrounding media satisfies the condition 
e'^/{ehvp) <C 1. Therefore, we claim that the energy difference between the first and zeroth Landau levels is always 
greater than the characteristic Coulomb attraction between the electron and the hole in the single graphene layer 
at any i? if e 3> /{hvp) « 2. Thus, applying perturbation theory with respect to weak Coulomb electron-hole 
attraction in graphene embedded in the GaAs microcavity (e = 12.9) is more accurate than for graphene embedded 
in the Si02 microcavity (e = 4.5). However, the magnetoexcitons in graphene exist in high magnetic field. Therefore, 
we restrict ourselves by consideration of high magnetic fields. 

Polaritons are linear superpositions of excitons and photons. In high magnetic fields, when magnetoexcitons may 
exist, the polaritons become linear superpositions of magnetoexcitons and photons. Let us define the superpositions 
of magnetoexcitons and photons as magnetopolaritons. It is obvious that magnetopolaritons in graphene are two- 
dimensional, since graphene is a 2D structure. The Hamiltonian of magnetopolaritons in the strong magnetic field 
is given by Eq. It can be shown that the interaction between two direct 2D magnetoexcitons in graphene with 
the electron on the Landau level 1 and the hole on the Landau level can be neglected in a strong magnetic field, in 
analogy to what is described in Ref. [17| for 2D magnetoexcitons in a quantum well. Thus, the Hamiltonian Htot © 
does not include the term corresponding to the interaction between two direct magnetoexcitons in a single graphene 
layer. So in high magnetic field there is the BEC of the ideal magnetoexcitonic gas in graphene. Therefore, in a single 
graphene layer is high mahnetic field we assume Hint = in Eq. ([3]). 

The binding energy f^**^ and effective magnetic mass ms of a magnetoexciton in graphene obtained using the first 
order perturbation respect to the electron-hole Coulomb attraction similarly to the case of a single quantum well 
are given by 

V 2 ers \/iTe^rB 

We obtain the effective Hamiltonian of polaritons by applying the standard procedure [i^ H^, [13, [Slj , when we 
diagonalize the Hamiltonian Htot © by using Bogoliubov transformations. If we measure the energy relative to the 
P = lower magnetopolariton energy {c / njtnr L'^^ — l^^^fll, we obtain the resulting effective Hamiltonian for trapped 

magnetopolaritons in graphene in a magnetic field. At small momenta a ^ 1 (Lp — hTic/n (ji!band — ^s^^ ) Sind 
weak confinement /3 ^ 1, this effective Hamiltonian is 

^- = E^^]^£^ + ^^M)^i>^p' (10) 



p 



and the effective magnetic mass of a magnetopolariton is given by 



(11) 



According to Eqs. pT|) and Q, the effective magnetopolariton mass Mcs increases with the increment of the magnetic 
field as i?^/^. Let us emphasize that the resulting effective Hamiltonian for magnetopolaritons in graphene in a 
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magnetic field for the parabolic trap is given by Eq. (jlOp for both physical realizations of confinement represented by 
case 1 and case 2. 

Neglecting anharmonic terms for the magnetoexciton-photon coupling, the Rabi splitting constant Qn can be 
estimated quasiclassically as 



1 /2 

\hnR\' ^ {f K^nt , EpM=(^^^^ , H^nt^-^^■A='^&■EpHO, (12) 



where a — {a^TO^y), and ay are Pauli matrices, Hint is the Hamiltonian of the electron-photon interaction cor- 
responding to the electron in graphene described by Dirac dispersion, Epho is the electric field corresponding to a 
single cavity photon, W is the volume of microcavity, uj is the photon frequency. The initial |i) electron state cor- 
responds to the completely filled Landau level and completely empty Landau level 1. The final |/) electron state 
corresponds to creation of one magnetoexciton with the electron on the Landau level 1 and the hole on the Landau 
level 0. The transition dipolc moment corresponding to the process of creation of this magnetoexciton is given by 
di2 = ers/^- Let us note that in Eq. the energy of photon absorbed at the creation of the magnetoexciton is given 
by hiij = El — So — \/2TLVF/'rB (we assume that 8^^^ ^ £i — Eq)- Substituting the photon energy and the transition 
dipole moment from into Eq. (|12p. we obtain the Rabi splitting corresponding to the creation of a magnetoexciton with 

1 /2 

the electron on the Landau level 1 and the hole on the Landau level in graphene: MIr = 2e {irhvprB / {V^fW)^ 

Thus, the Rabi splitting in graphene is related to the creation of the magnetoexciton, which decreases when the 
magnetic field increases and is proportional to Therefore, the Rabi splitting in graphene can be controlled by 

the external magnetic field. It is easy to show that the Rabi splitting related to the creation of the magnetoexciton, the 
electron on the Landau level and the hole on the Landau level — 1 will be exactly the same as for the magnetoexciton 
with the electron on the Landau level 1 and the hole on the Landau level 0. Let us mention that dipole optical 
transitions from the Landau level —1 to the Landau level 0, as well as from the Landau level to the Landau level 1, 
are allowed by the selection rules for optical transitions in single- layer graphene [53] . 

Although Bose-Einstein condensation cannot take place in a 2D homogeneous ideal gas at non-zero temperature, as 
discussed in Ref. [37j . in a harmonic trap the BEC can occur in two dimensions below a critical temperature T^. In a 
harmonic trap at a temperature T below a critical temperature (T < T^), the number No{T, B) of non- interacting 
magnetopolaritons in the condensate is given in Ref. (STj . Applying the condition = 0, and assuming that the 

magnetopolariton effective mass is given by Eq. (|lip . we obtain the BEC critical temperature ri"' for the ideal gas 
of magnetopolaritons in a single graphene layer in a magnetic field: 



where N is the total number of magnetopolaritons, g^^'^^^ and gi'^-'''''^ are the spin and graphene valley degeneracies 
for an electron and a hole, respectively, ks is the Boltzmann constant. At temperatures above Ti^\ the BEC of 
magnetopolaritons in a single graphene layer does not exist. T^^'^ /\fN as a function of magnetic field B and spring 
constant 7 is presented in Fig. [21 In our calculations, we used glf' — gi'^'' — gi'^'' — gi'^^ — 2. According to Eq. (|13p . 
the BEC critical temperature ri"-* decreases with the magnetic field as B~^/^ and increases with the spring constant 
as 7^/^. Note that we assume that the quality of the cavity is sufficiently high, so that the time of the relaxation to 
the Bose condensate quasiequilibrium state is smaller than the life time of the photons in the cavity. 

Above we discussed the BEC of the magnetopolaritons in a single graphene layer placed within a strong magnetic 
field. What would happen in a multilayer graphene system in a high magnetic field? Let us mention that the 
magnetopolaritons formed by the microcavity photons and the indirect excitons with the spatially separated electrons 
and holes in different parallel graphene layers embedded in a semiconductor microcavity can exist only at very low 
temperatures fc^T <C hrtn. For the case of the spatially separated electrons and holes, the Rabi splitting fi/j is very 
small in comparison to the case of electrons and holes placed in a single graphene layer. This is because Qr ~ di2 
and the matrix element of magnetoexciton generation transition di2 is proportional to the overlapping integral of the 
electron and hole wavefunctions, which is very small if the electrons and holes are placed in different graphene layers. 
Therefore, we cannot predict the effect of relatively high BEC critical temperature for the electrons and holes placed 
in different graphene layers. 
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FIG. 2: The ratio of the BEG critical temperature to the square root of the total number of magnetopolaritons Tc^^ /\/N as 
a function of magnetic field B at different spring constants 7. We assume the environment around graphene is GaAs with 
e = 12.9. 



IV. SUPERFLUIDITY OF MAGNETOEXCITONS IN BILAYER GRAPHENE 

We consider two parallel graphene layers (GLs) separated by an insulating slab of dielectric (for example, Si02) [s^ ]. 
The spatial separation of electrons and holes in different GLs can be achieved by applying an external electric field. 
Besides, the spatially separated electrons and holes can be created by varying the chemical potential by applying a 
bias voltage between two GLs or between two gates located near the corresponding graphene sheets. The equilibrium 
system of local pairs of spatially separated electrons and holes can be created by varying the chemical potential by 
using a bias voltage between two GLs or between two gates located near the corresponding graphene sheets (case 1) (for 
simplicity, we also call these equilibrium local e-h pairs in magnetic field B as indirect magnetoexcitons) . In the case 
1 a magnetoexciton is formed by an electron on the Landau level 1 and hole on the Landau level —1. Magnetoexcitons 
with spatially separated electrons and holes can be created also by laser pumping (far infrared in graphene) (case 2) 
and by applying perpendicular electric field as for CQWs 24, 25, 27]. In the case 2 a magnetoexciton is formed by an 
electron on the Landau level 1 and hole on the Landau level 0. We assume the system is in quasiequilibrium state. 
Below we assume the low-density regime for magnetoexcitons, i.e. magnetoexciton radius a < n~^/^, where n is the 
2D magnetoexciton density. 

In a strong magnetic field at low densities, n <C (r^ — [he/ {eB))^^'^ is the magnetic length, — e is the electron 
charge, c is the speed of light), indirect magnetoexcitons repel as parallel dipoles, and we have for the pair interaction 
potential: 

where D is the interlayer separation, e is the dielectric constant of the insulator between two layers, R.i(2) are radius 
vectors of the center of mass of two magnetoexcitons. Since typically, the value of r is tb, and P <C h/rs in this 
approximation, the effective Hamiltonian Hmex in the magnetic momentum representation P in the subspace the lowest 
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Landau level has the same form (compare with Ref. fl3l) as for two-dimensional boson system without a magnetic 
field, but with the magnetoexciton magnetic mass ttib (which depends on B and D; see below) instead of the exciton 
mass {M = + mh)^ magnetic momenta instead of ordinary momenta. We can obtain the effective Hamiltonian 
for bilayer graphene without a confinement if we keep considering only two graphene layers in the magnetoexciton 
effective Hamiltonian without a trap ([3]) : 

ffmea; = X!'^o(-P)apap + ^ ^ <^Pl , P2 P3 , P4 ^ Op^ OPa > (15) 

P Pi.P2,P3,P4 

where the matrix element (Pi,P2 



t/ P3,P4^ is the Fourier transform of the pair interaction potential U{R) — 

e^D^ /eB? and for the lowest Landau level we denote the spectrum of the single exciton eo{P) = £oo(P)- For an isolated 
magnetoexciton on the lowest Landau level at the small magnetic momenta under consideration, eo(P) ~ ^'^/(Stib), 
where is the effective magnetic mass of a magnetoexciton in the lowest Landau level and is a function of the 
distance D between e - and h - layers and magnetic field B (see Ref. [2H). In strong magnetic fields at ^ the 
exciton magnetic mass is mB{D) = eD^/{e^r%) for the QWs [2l| and mB{D) = eD^ /{ie'^rj^) for the GLs ^52j . 

We study the magnetoexciton-magnetoexciton scattering applying the theory of weakly-interacting 2D Bose-gas [H, 
[Tij . The chemical potential fi of two-dimensional dipole magnetoexcitons in graphene bilayer system, in the ladder 
approximation, has the form (compare to Refs. [igI. [lOj ) : 

^ ^ TT^^n 

where s = 8 is the spin and valley degeneracy factor for a magnetoexciton in graphene bilayer, n is the 2D density of 
magnetoexcitons. 

At small momenta the collective spectrum of magnetoexciton system is the sound-like e{p) = CgP (cj = -y/ /i/(2tob) 
is the sound velocity) and satisfied to Landau criterium for superfluidity. The density of the superfluid component 
ns{T) for two-dimensional system with the sound spectrum can be estimated as:[3^ 

, 3C(3) k%T^ , , 

where C,{z) is the Riemann zeta function, and ^(3) ~ 1.202. The second term in Eq. (fT7|) is the temperature dependent 
normal density taking into account gas of phonons ("bogolons") with dispersion law e(p) = ^ fi/ {2mB)p, H is given 
by Eq.lIlD. 

In a 2D system, superfluidity of magnetoexcitons appears below the Kosterlitz-Thouless transition temperature [sj : 
Tc = TTns{T)/{2mB)- The dependence of Tc on B and D for the cases 1 and 2 is represented in Fig. [3] (since in the 
case 1 the binding energy two times higher and the effective magnetic mass is two times smaller than in the case 
2 [s^l , the magnetoexcitons in the case 1 are expected to be twice more stable and Tc in the case 1 is expected to be 
approximately twice higher than in the case 2 at fixed n, D and B). The temperature Tc for the onset of superfluidity 
due to the Kosterlitz-Thouless transition at a fixed magnetoexciton density decreases as a function of magnetic field 
B and interlayer separation D. This is due to the increased tub as a functions of B and D. The Tc decreases as 
B'^/"^ at £) < rs or as B'"^ when D > r^. 



V. INSTABILITY OF DIPOLE MAGNETOEXCITONS AND SUPERFLUIDITY OF 
MAGNETOBIEXCITONS IN GRAPHENE SUPERLATTICES 

We consider magnetoexcitons in the superlattices with alternating electronic and hole GLs. We suppose that 
recombination times can be much greater than relaxation times Tr due to small overlapping of spatially separation 
of electron and hole wave functions in GLs. In this case electrons and holes are characterized by differe nt q uasi- 
equilibrium chemical potentials. Then in the system of indirect excitons in superlattices, as in CQW [H [U, the 
quasiequilibrium phases appear. No external field applied to a slab of superlattice is assumed. If "electron" and 
"hole" quantum wells alternate, there are excitons with parallel dipole moments in one pair of wells, but dipole 
moments of excitons in another neighboring pairs of neighboring wells have opposite direction. This fact leads to 
essential distinction of properties of e — /i system in superlattices from one for coupled quantum wells with spatially 
separated electrons and holes, where indirect exciton system is stable due to dipole-dipole repelling of all excitons. 
This difference manifests itself already beginning from three-layer e — h — e or h~e~h system. We assume that 
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FIG. 3: Dependence of Kosterlitz-Thouless transition temperature Tc — Tc{B) (in units K) versus magnetic field for two-layer 
graphene separated by Si02. with ej, = 4.5. The magnetoexciton density n — 4 x 10^^ cm~^. Different interlayer separations 
D are chosen for the case 1: D — 30nm (solid curve), D = 28nm (dotted curve), D = 27nm (dashed curve). For the case 2: 
D — 30nm (dashed-dotted curve), D = 28nm (thin solid curve), D = 27nm (thin dotted curve). 

alternating e — h — e layers can be formed by independent gating with the corresponding potentials which shift chemical 
potentials in neighboring layers up and down or by alternating doping (by donors and acceptors, respectively). 

Let us show that the low-density system of weakly interacting two-dimensional indirect magnetoexcitons in super- 
lattices is instable, contrary to the two-layer system in the CQW. At the small densities n the system of indirect 
excitons at low temperatures is the two-dimensional weakly nonideal Bose gas with normal to wells dipole moments 
d in the ground state {d = eD, D is the interlayer separation). In contrast to ordinary excitons, for the low-density 
spatially indirect magnetoexciton system the main contribution to the energy is originated from the dipole-dipole 
interactions U- and U+ of magnetoexcitons with opposite (see Fig. [1]) and parallel dipoles, respectively. 

The behavior of the potential energies U+ {R) and C/_ (i?) as the functions of the distance between two excitons R 
is shown in Fig. U) We suppose that D/R ^ 1 and L/R 1, where L is the mean distance between dipoles normal 
to the wells. We consider the case, when the number of quantum wells k in superlattice is restricted k ^ l/{Dy/Trn). 
This is valid for small k or for sufficiently low exciton density. 

We can obtain the effective Hamiltonian for graphene superlattice without a confinement if we keep considering only 
the magnetoexciton effective Hamiltonian ^ without a trap. Let us apply the Bogolubov approximation to analyze 
a stability of the ground state of the weakly nonideal Bose gas of indirect excitons in superlattices. We assume U+ 
and U- are the 2D Fourier images of U+{R) and U-{R) at P = 0, respectively, and S is the surface of the system. 
Let us mention that the appropriate cut-off parameter for this Fourier transform is the classical turning point of the 
dipole-dipole interaction. Let us mention that the appropriate cut-off parameter for this Fourier transform is the 
classical turning point of the dipole-dipole interaction. Note that the cut-off parameter Rq for the potential U+{R) is 
much greater than for U-{R) (the cut-off parameters i?o for the both potentials can be represented in Fig. |4]by the 
points where the curves corresponding to U+{R) and U^{R) are crossed by the chemical potential /x^ represented by 
the horizontal straight line placed right above but close to L/_|_(_)(i?) — 0). Therefore, we claim that U+ > 0, [/_ < 0, 
and \U-\ > \U+\. 

Applying the unitary Bogoliubov transformations to the magnetoexciton operators Op, Op, bp, and bp, we diago- 
nalize the Hamiltonian Htot in the Bogoliubov approximation 55]. Finally, we obtain 

Htot = eb)("pap + fipPp) (18) 

with the spectrum of quasiparticles e{P): 

el(P) - el{P) + 2nU+e^{P), 

4iP) = el{P) + 2n{U+ + U^)eo{P) . (19) 

Since U+ > and /7_ < 0, we have el{P) > S2{P) at P > 0. Therefore, at low temperatures the quasiparticles 
only with the spectrum S2{P) will be excited, since the excitation of these quasiparticles requires less energy than 
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FIG. 4: The potential energy of the interaction of indirect magnetoexcitons with parallel U+{R) and opposite U-{R) dipoles, 
located in neighboring pairs of GLs (in units of the binding energy of the indirect magnetoexciton E^x = /eD), as a function 
of the distances R between magnetoexcitons along the GLs (in units of D). 

for the quasiparticles with the spectrum el{P). Since [/+ + U- < 0, it is easy to see from Eq. (IT9|) that for the 
small momenta P < Am.Bn\U+ + U-\ the spectrum of excitations becomes imaginary. Hence, the system of weakly 
interacting indirect magnetoexcitons in the slab of the superlattice is unstable. It can be seen that the condition of the 
instability of magnetoexcitons becomes stronger as magnetic field higher, because tub increases with the increment 
of the magnetic field, and, therefore, the region of P resulting in the imaginary collective spectrum increases as B 
increases. 

As the ground state of the system we consider the low-density weakly nonideal gas of two-dimensional indirect 
magnetobiexcitons, created by indirect magnetoexcitons with opposite dipoles in neighboring pairs of wells (Fig. [1]). 
The mean dipole moment of indirect magnetobiexciton is equal to zero. However, the quadrupole moment is nonzero 
and equal to Q = 3eZ3^ (the large axis of the quadrupole is normal to quantum wells/graphene layers). So indirect 
magnetobiexcitons interact at long distances D as parallel quadrupoles: U{R) — 9e^Z3^/(ei?^). 

We apply the theory of weakly-interacting 2D Bose gas [lB| to study the magnetobiexciton-magnetobiexciton 
repulsion. The chemical potential /i of two-dimensional biexcitons, repulsed by the quadrupole law, in the ladder 
approximation, has the form (compare to Refs. (l6l. [l9j): 

^ m^l0g[n4/3g2/3/(8^(18^^e2254)2/3„^^^)]- ^ > 

where Uhex — n/St is the density of magnetobiexcitons in graphene layers and we considered that the magnetic mass 
of magnetobiexciton is twice of the magnetic mass of magnetoexciton, i.e. 2mB- 

At small momenta the collective spectrum of magnetobiexciton system is the sound-like e{p) — CsP (c^ — ^ [ij {2mB) 
is the sound velocity) and satisfied to Landau criterion for su per fluidity. The density of the superfluid component 
ns{T) for two-dimensional system with the sound spectrum is [32] 

n5(r) = n,..-^^. (21) 
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FIG. 5; Dependence of the Kosterlitz-Thouless transition temperature Tc — Tc{B) for the superlattice consisting of the QWs 
for GaAs/AlGaAs, e = 13; and for GLs separated by the layer of Si02 with e = 4.5 on the magnetoexciton density n at D = 10 
nm at different magnetic fields. The solid, dashed and thin solid curves for the QWs, dotted, dashed-dotted and thin dotted 
curves for the GL at B: B = 20 T, B = 15 T and B = 10 T, respectively. 

In a 2D system, superfluidity of magnetobiexcitons appears below the Kosterlitz-Thouless transition temperature 
Tc = TTns(T)/{AmB)- Employing ns(T) for the superfluid component, we obtain an equation for the Kosterlitz- 
Thouless transition temperature Tc Tc = Tms{T) / {Atub) ■ The dependence of Tc on the density of magnetoexcitons 
at different magnetic field B for superlattice consisting of quantum wells and graphene layers is represented on Fig. [S] 
Let us mention, that we have apply the same effective Hamiltonian to describe magnetobiexcitons in high magnetic 
field in the superlattice of QWs and GLs with the only difference in the effective magnetic mass of magnetoexciton 

771b. 

According to Fig. [5l the temperature Tc for the onset of superfluidity due to the Kosterlitz-Thouless transition at a 
fixed magnetoexciton density decreases as a function of magnetic field B and interlayer separation D. This is due to 
the increased effective magnetic mass 777,5 of magnetoexcitons as a functions of B and D. The Tc decreases as B~^^^ 
at D <C or as B~^ when D tb- According to Fig. [51 the Kosterlitz-Thouless temperature Tc is higher for the 
superlattice consisting of graphene layers than for the superlattice consisting of the quantum wells, and this difference 
is as stronger as the magnetic field is smaller. 

VI. CONCLUSIONS 

We have obtained the effective Hamiltonian of the quasiparticles in graphene structures in high magnetic field: 
indirect magnetoexcitons in graphene bilayer, magnetobiexcitons in graphene superlattices and magnetopolaritons in a 
graphene layer embedded in optical microcavity. It was shown that the gas of magnetoexcitons in graphene superlattice 
is instable due to the attraction between magnetoexcitons with parallel dipoles, while the system of magnetobiexcitons 
in the graphene superlattice is stable. We have shown that, while the quasiparticles in GLs and QWs in high magnetic 
field are described by the same effective Hamiltonian with the only difference in the effective magnetic mass of 
magnetoexciton. This difference is caused by the four-component spinor structure of magnetoexciton wave function 
in GL, while magnetoexciton in a QW and CQWs is characterized by the one-component scalar wave function. Besides, 
we show that the magnetoexciton system graphene bilayer and magnetobiexciton system in graphene superlattice can 
be described as a 2D weakly-interacting Bose gas, which is superfluid below the Kosterlitz-Thouless phase transition 
temperature. We have calculated the density of the superfluid component and the Kosterlitz-Thouless temperature for 
the systems of magnetoexcitons and magnetobiexcitons as functions of magnetic field B, interlayer separation D, and 
magnetoexciton density n. In contrast to the magnetoexcitons in graphene bilayer and magnetobiexcitons in graphene 
superlattice, the magnetopolaritons in GL embedded in an optical cavity in the limit of high magnetic field is ideal 
Bose gas without interparticle interactions, and, therefore, magnetopolariton gas is not superfluid. However, there is 
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is BEC at the temperatures below critical one in this system in a trap. We have calculated the critical temperature 
of magnetopolariton BEC in GL embedded in an optical microcavity in a trap as a function of magnetic field and the 
curvature of the trap 7. Note that taking into account the virtual transitions of electrons and holes between Landau 
levels results in weak (at large e) interactions between magnetoexcitons [ij'l. In turn, this leads to the possibility of 
the superfluidity of the magnetopolariton system. 
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